Rate of photon production from hot hadronic matter 
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Thermal photon emission rates from hot hadronic matter are studied to order e g , where g indi- 
cates a strong-interaction coupling constant. Radiative decay of mesons, Compton and annihilation 
processes for hadrons, and bremsstrahlung reactions are all considered. Compared to the standard 
rates from the literature, one finds two orders of magnitude increase for low photon energies stem- 
ming mainly from bremsstrahlung and then a modest increase (factor of 2) for intermediate and high 
energy photons owing to radiative decays for a variety of mesons and from other reactions involving 
strangeness. These results could have important consequences for electromagnetic radiation studies 
at RHIC. 



PACS numbers: 25.75.-q, 12.38.Mh 

Photons carry undistorted detail on strong-interaction 
dynamics since their mean free paths in hot hadronic 
matter are much larger than typical system sizes. For 
this reason, photon measurements in high energy nuclear 
collisions are deemed the best thermometers of the high- 
est energy densities reached thereP,|3l- While pions and 
protons scatter several times before leaving the collision 
zones, photons do not, and therefore originate from the 
entire spacetime hadronic source, and in their asymp- 
totic momentum configurations. Quantitative studies of 
photon mean free paths in hadronic matter have indeed 
revealed A ~ 2x10^ fm for a temperature 200 MeVy. 

Since production rates are increasing functions of tem- 
perature, the inverse slope on a transverse momentum 
spectrum, for instance, is a good measure of the kinetic 
properties of the early (thermal) hadronic source. And 
yet, photons are produced during the entire reaction. 
Pre-equilibrium dynamics contribute, parton dynamics 
contribute, and then, hadron dynamics might be con- 
sidered as a background to the more sought-after quark 
gluon plasma contributions 0, 0|. For quantitative in- 
terpretation of photon measurements at the Relativistic 
Heavy Ion Collider (RHIC), a complete determination of 
all components, including the hadronic background, is 
clearly a prerequisite. 

The early comparison of photon production from quark 
gluon plasma versus hadronic matter by Kapusta et al.^ 
revealed that the two phases have equal production rates 
at fixed temperature. Technically, the advancements 
here were made possible through an application of hard- 
thermal-loop resummation methods 0, and effective 
hadronic field theory. This physically appealing result 
is however, somewhat discouraging in terms of allowing 
distinction between the two phases using photons. The 
equal-luminous property of the two phases has survived 
even after next-to-leading-order (NLO) QCD processes 
have been considered I fll . NNLO parton processes have 
even been analvzed |ld IllL H^. and the new rate is still 
consistent with the hadron rates. The inevitable conclu- 



sion is that photon production as a QGP diagnostic is not 
as definitive as was originally hoped. However, with the 
tremendous recent advancements in the determination of 
QCD processes — to three-loop order, including annihila- 
tion with scattering and also bremsstrahlung with multi- 
ple scattering interferences fl^ . it behooves theory to 
put the hadronic formalism on a more equal footing. The 
purpose of this paper is to report on photon production 
from hadronic matter computed at NNLO in the strong- 
interaction coupling constant, namely, to order e^g*. 

The photon production rate to lowest order in the elec- 
tromagnetic coupling and to all orders in the strong cou- 
pling can be compactly written as[l5L llfi[ll7l | 
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where the retarded photon self-energy includes all possi- 
ble hadron topologies — that is, to arbitrary loop order. 
The imaginary part corresponds to emission and absorp- 
tion, of which the former is relevant here. Absorption 
was studied in Ref. and foimd to be negligible. At the 
one-loop level, the imaginary part corresponds to such 
lowest-order reactions as w ^ irj and p ^ nj. There 
is of course a long list of other resonance decays to be 
considered, as discussed below. Typically, in order to 
quantify production rates, a model for the strong inter- 
action (e.g. an effective Lagrangian) is necessary. How- 
ever, the one-loop reactions have the desirable feature 
that the rate can be computed in a model independent 
way, requiring only the measured radiative decay rates. 
If such measurements are not available, then vector me- 
son dominance plus a chiral Lagrangian is an appropriate 
approach. 

For the process a — s- 1 + 7, where a and 1 are mesons 
with allowed quantum numbers and masses, the thermal 
rate can be written as 
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TABLE I: Empirical radiative decay widths for low-lying res- 
onances. 



Process 



Width (keV) 



p(770) 
cj(783) - 
K*{892) 
<^(1020) ^ -q-y 
61 (1235) ^ -K-y 
01 (1260) ny 
a2(1320) -Kj 
A'i(1270) ^ K-/ 



117.9 
734.2 
83.5 
55.3 
227.2 
640±240 
286 
70" 



"Vector meson dominance was used to estimate this value. 



X [1 + fBEiEa-E^)], 



(2) 



where M is the spin and isospin degeneracy for the par- 
ent, Eq is the photon energy in the rest frame of the 



parent, E„ 



maiE^ + E^)/2E^Eo, and fsE is the 



Bose-Einstein distribution to account for quantum statis- 
tics and medium-enhancement for species 1. If baryons 
were studied, e.g. A — > A^7, then Fermi-Dirac distribu- 
tions would be appropriate and a minus sign in the sec- 
ond factor of the integrand would suppress phase-space 
occupation for species 1. The above integral is typically 
reported in the literature, and then a numerical exercise 
ensues. However, there is no need to leave it here since 
the integral has a closed-form expression with the result 
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Empirical radiative decay widths for several low-lying 
resonances are available! IS] and included in Table The 
resulting thermal decay rates are displayed in Fig. ^ 

Next-to-leading-order contributions enter at order e^g^ 
and come from two-loop self-energy structures. Cou- 
pling strength and phase space arguments suggest that 
TT/j — > 7r7, TTTT — > p7 and p tttt^ will be important con- 
tributors. The first channel, which is the Compton type 
process, will be considered here while the second and 
third will be discussed later in the higher orders, since 
they are essentially the on-shell components of the full, 
higher order contributions. For the Compton process, a 
particular interaction Lagrangian was assumed in Ref. 
describing ir-p dynamics, and then calibrated to p — )■ tttt 
decay. In terms of practical kinetic theory formalism, the 
rate for the general reaction a + 6^1+7is 

^<Pp^ J {2TTY2Ea {2TTf2Eb {2TTf2El^^ 

\M\''{2T:f6^{pa+Pb-pi-p^)i^) 
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FIG. 1: Photon production rates from lowest-order reactions 
of radiative meson decays at fixed temperature T = 200 MeV. 



A numerical task remains: That is, to integrate Eq. © 
over the appropriate phase space. The result for Trp — > 7r7 
is shown in Fig. [5] 

There are many other contributions at the same order 
with roughly the same phase space properties. Certain 
strangeness reactions turn out to be strong contributors, 
as can be seen from the results plotted in Fig. |21 The 
coupling constants are fitted to empirical hadronic decays 
where available, and estimated using an SU(3) symmetry 
in a chiral Lagrangian in the absence of measurements. 
The strongest contributor involving strangeness seems to 
be TiK* Kj, although there are certainly others of 
some significance. Strangeness reactions of this type go 
beyond what has been done before 20] . Also, these two- 
loop reactions and others are currently being studied in a 
parallel investigation, broadened in scope to explore the 
issue of form factors 21]. 

A remark is in order here regarding the role of the Oi 
meson. In the reaction irp — > 7r7, the ai appears as an in- 
termediate state, both in the s channel and the u channel. 
Early estimates reported that the s channel dominated 
all others, even dominated tt exchangel^^. Later stud- 
ies using a hidden local symmetry in a chiral Lagrangian 
were less clear regarding the role of the ai since no unique 
rate could be establishedjl^. However, recent improve- 
ments describing ai dynamics with particular concern for 
the D/S ratio in the oi — + 7rp scattering amplitude show 
a diminished iinportance of the role of the ai for pho- 
ton production|2J]. In Fig. ^ the oi radiative decay has 
been included while in Fig. [21 only the 7r-exchange result 
for the Compton process np — > 7r7 is shown. This has 
the slight incomplete feature of ignoring the u channel 



3 



10 
10" 
\ 10" 
k 10" 
qT-IO" 

n 
"O 

^ 10" 



10 r 



10 




T = 200 (MeV) 
I I I L 







1 2 
(GeV) 



FIG. 2: Thermal photon production rate from important 
channels at order e^g^. 
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FIG. 3: Bremsstrahlung contributions (order e^g"^) to thermal 
photon production rate at T = 200 MeV. 



oi exchange and coherence effects, but at the same time 
has the advantage of removing any possibility for double 
counting. 

Motivated partly by the realization that the QCD pro- 
cesses and the resulting photon rates have been computed 
with bremsstrahlung included (and the hadron rates up 
to now have not), one moves next to three- loop topolo- 
gies. Three-loop self-energy graphs involve four on-shell 
hadron lines with a variety of possible kinematics, plus 
the photon. Some of these include processes of the type 
a + h — > I-I-2-I-7: hadronic elastic scattering with 
bremsstrahlung. Based on general principles alone, these 
contributions are expected to be strong for low photon 
energies and quickly become less important for increas- 
ing photon energy. Indeed, they will even diverge as 
^ Q. When bremsstrahlung processes are included, 
there is immediately a danger of double counting. For 
example, tttt — > p7 and p tttt"/ are already part of the 
full TTTT — > 7r7r7 coherent calculation with an intermedi- 
ate p meson. Even though the two separate reactions 
were generated from two-loop self-energy topologies (re- 
stricted, however, to on-shell vector mesons), they must 
not be added on top of the bremsstrahlung contributions. 
That would be double counting. The strategy taken here 
is the following. When two-loop contributions are sub- 
sets of three-loop contributions, the separate two-loop 
contributions will not be included. 

A treatment of the bremsstrahlung mechanisms in 
terms of complete phase space is a bit tedious, and so 
to first assess their relative importance, a soft-photon 



approximation is adopted here. The rate is thenj25l| 

dz A(z^r^, TO^, ml) Ki{z) 



E. 



da 



R2{s,ml,ml) ' 



(5) 



where 



nT — E^ 



M2 + E^, with M = max(r 



nib, mi + m2), where the full squared amplitude is ap- 
proximated as \M\'^ « \Mo\'^i—J'^), where Mo is the 
elastic scattering amplitude, J'^ is the hadronic current 
including the fundamental charge e, Ki is the modified 
Bessel function of order 1, and finally, where the ratio of 
the two-body phase space for S2 = s — 2y/sE^ and s re- 
spectively, is a phase-space correction factor to minimize 
error at finite E-y. 

Bremsstrahlung reactions are shown in Fig. |31 for pi- 
ons, rho mesons, K and K* mesons. Interactions are 
once again modeled with an SU(3) chiral Lagrangian 
with coupling constants fitted by respecting hadronic 
phenomenology. As expected, the partial rates are very 
strong at low photon energy. But the striking feature 
is the strength surviving even for intermediate photon 
energies. 

Two cautionary remarks must however be made. First, 
Eq. (O does not include multiple scattering interference 
effects of Landau-Pomeranchuk-Migdal, which will re- 
duce the rates when the photon energy is less than the 
inverse of the mean time between strong interactions and 
second, the soft photon approximation becomes less re- 
liable with increasing photon energy. Still, one does not 
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gies reported here in the new rate as compared with the 
old rate is not quantitatively determined. However, the 
general conclusions are expected to be robust. It will 
therefore become extremely important to assess such ef- 
fects in future studies, as well, to incorporate these new 
rates into dynamical or statistical models for the pur- 
poses of comparison with experimental results at RHIC 
and interpretation regarding QGP diagnostics. 
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FIG. 4: Comparison of conventional hadronic rate to or- 
der e^g^ ("old rate" as parameterized by Nadeau et al. |26ll') 
with an updated rate including additional radiative de- 
cays, strangeness reactions, and bremsstrahlung ( "new rate" ) 
which goes to order e^g'^. The old rate includes tttt p7, 
p 7r7r7, 7rp — > 7r7 and lu — > Tvy. 



expect more than a factor of 2-4 suppression for ^ 
0.5 GeV. 

Finally, the most useful comparison one can make re- 
garding the relative importance of the higher-order con- 
tributions is to add up all the partial rates from radiative 
decays (order e^), Compton processes (order e^g^), and 
from bremsstrahlung (order e^g'*) and then compare to 
the conventional hadronic rate from the literature (which 
includes uj — s- n'^j, tttt — s- p7, p — s- mrj and irp — > irj). 
This comparison is reported in Fig. ^ A factor of 2 in- 
crease is seen over most of the intermediate and high 
energy region up to 3 GeV. Below 1 GeV photon energy, 
the increase is rather large, approaching two orders of 
magnitude around 200 MeV photon energy. 

One might summarize this study as follows. Pho- 
ton production rates in hot hadronic matter have been 
investigated up to next-to-next-to-leading order in the 
strong coupling. This means contributions to order e^g^ 
have been included, most notably, elastic scattering of 
mesons with accompanying bremsstrahlung. In addi- 
tion, strangeness reactions were included and found to 
be somewhat important. Quantitatively, the new photon 
rate with strangeness and with NNLO hadronic processes 
is two orders of magnitude greater at low photon energies 
and a factor of approximately 2 greater at intermediate 
and high energies as compared with the standard photon 
production rate from the literature [2^. 

Since multiple scattering interferences have not yet 
been included, the precise increase at low photon ener- 
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